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gravity with a dynamical preferred frame 
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A generally covariant extension of general relativity (GR) in which a dynamical unit timelike vec- 
tor field is coupled to the metric is studied in the asymptotic weak field limit of spherically symmetric 
static solutions. The two post-Newtonian parameters known as the Eddington-Robertson-Schiff pa- 
rameters are found to be identical to those in the case of pure GR, except for some non-generic 
values of the coefficients in the Lagrangian. 
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I. INTRODUCTION 

Over the past several years quantum gravity consider- 
ations have led a number of researchers to contemplate 
violations of Lorentz invariance. Such violations are usu- 
ally studied within the realm of particle physics, but they 
would also have implications for gravitation. Fixed back- 
ground tensors breaking Lorentz symmetry also break 
the general covariance of general relativity (GR), gener- 
ically resulting in over-determined equations of motion. 
Together with a theoretical bias for preserving general 
covariance, this leads one to consider Lorentz breaking 
introduced by dynamical tensor fields. 

The background value of a given dynamical tensor 
field "spontaneously" breaks Lorentz symmetry, either 
because of a potential that has a Lorentz violating (LV) 
minimum, or because the tensor is somehow constrained 
not to vanish. The simplest situations arise from scalar 
fields whose non-zero gradient is an LV vector, and vector 
fields. Vector fields without any potential or constraint 
were considered in the early f970's by Will, Nordvedt, 
and Hellings in the spirit of alternate theories of grav- 
ity 1^ IM S' Vector fields with a potential leading 
to LV were studied by Kostelecky and Samuel 8] mo- 
tivated by string theory considerations, and by Clayton 
and Moffat 1^ motivated by the notion that a dynami- 
cally varying speed of light might solve some cosmological 
problems (see also IEH)- 

A vector field constrained to have a fixed timelike or 
spacelike length breaks local Lorentz symmetry in every 
configuration, much as a nonlinear sigma model sponta- 
neously breaks gauge invariance. In the timclikc case the 
residual symmetry is the 3D rotation group, while in the 
spacelike case it is the 2+1 Lorentz group. A particu- 
larly simple example of such a theory was considered by 
Kostelecky and Samuel and also studied by Jacobson 
and Mattingly 0. In this example, the action for the 
covariant vector field is just the square of the exterior 
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derivative. Like the Maxwell action for the vector poten- 
tial, this is independent of the connection components. 

The most general action with up to two derivatives 
of a unit vector field has four terms. Three of these 
terms were included in the original study of Will and 
Nordvedt P], (which did not include any constraint on 
the length of the vector) and all four were written down 
and studied using the tetrad formalism by Gasperini , 
who broke the local Lorentz symmetry by including in the 
action terms referring to a fixed "internal" unit timelike 
vector. The same theory in the metric formalism was 
written down in Q , and the linearized wave solutions in 
one special case (cor resp onding to the case focused on in 
1]) were reported in ICrl. The general constrained vector- 
tensor theory is quite complicated due to the derivative 
coupling to the vector field which includes connection 
components. Therefore the action for the vector in fact 
modifies the kinetic terms for the metric as well. Thanks 
to the unit timelike constraint, the vector has only three 
degrees of freedom, all corresponding to spacelike varia- 
tions. Thus, unlike in other vector theories without gauge 
invariance, problems with negative energy modes need 
not arise. 

We are interested in the observational signatures and 
constraints on the parameters in the general constrained 
vector-tensor theory. These can be obtained from the 
parametrized post-Newtonian (PPN) parameters, wave 
phenomena, and strong field effects. The PPN param- 
eters for the general unconstrained vector-tensor theory 
were found by Will i^]. Those results do not directly 
apply to the constrained case, hence the analysis must 
be carried out anew. In this paper we begin that pro- 
cess, restricting at first to the static PPN parameters, 
i.e. the Eddington-Robertson-Schiff parameters /? and 
7, which are the only ones that do not vanish when the 
isolated system is at rest with respect to the asymptotic 
rest frame defined by the vector field. We find that for 
generic values of the coefficients in the Lagrangian these 
parameters take the same values as in GR. This indicates 
that to observationally bound the coefficients in the la- 
grangian one must consider higher order PPN contribu- 
tions, preferred frame effects associated with the motion 
of the solar system relative to the asymptotic rest frame 
of u"-, and/or radiation or other effects. 

We use metric signature (H ) and units with c — 1. 



2 



II. ACTION AND FIELD EQUATIONS 

Taking the viewpoint of effective field theory, we con- 
sider the action as a derivative expansion, keeping all 
terms consistent with diffeomorphism symmetry. The 
most general Lagrangian scalar density involving the 
metric gab and preferred frame unit vector m° with two 
or fewer derivatives is £ = y/—gC, with 

C = -R- i^"^,„(Vaw'")(V,^/") - Xigabu'^u" - 1) (1) 

where 

K'^'mn = Cig''''gmn + C2S'^S'^+C3S!^,5>^ + CiU^u'>gran. (2) 

Stationarity under variation of the Lagrange multiplier A 
constrains the preferred frame vector u° to be unit time- 
like. We have omitted in the lagrangian terms that would 
vanish when the unit constraint is satisfied. The unit 
timelike vector u"^ is present everywhere in spacetime in 
every field configuration and specifies a locally preferred 
rest frame. It can be thought of as the four-velocity of a 
ubiquitous fluid and hence is naturally called the aether 
field. Wc sometimes use the term aether theory to refer 
to the theory described by the Lagrangian (Q). 

A general class of Lagrangians for vector-tensor theo- 
ries with an unconstrained vector was parametrized by 
Will and Nordvedt Their term uju°-UaR does not 

appear here because of the unit constraint, and our C4 
term does not appear there. The relation between the 
parameters in Ref. 4] (neglecting lo) and ours (neglect- 
ing C4) is 

T = -{ci + C2 + C3) ci = 2e - r 

ry = -C2 C2 = -?7 (3) 

e = -(c2 + C3)/2 C3 = r/-2e 

In order to agree with observations the dimensionless 
coefficients ci_2.3,4 must presumably be fairly small com- 
pared to unity. Special cases of this action are identified 
in Table |3 Note that the "Einstein-Maxwell-like" case 
is a sub-case of Hellings-Nordtvedt. This case has an 
extra gauge symmetry and was disfavored in Q on ac- 
count of the gradient singularities that generally develop 
in the vector field. Another notable sub-case of Hellings- 
Nordtvedt is ci = C4 = 0, C2 -I- C3 = 0, which is also one 
of the theories equivalent to GR via a field redefinition. 

The metric equation with no matter source (other that 
the aether field) can be written in the form 

Gab = Tab (4) 

where Gab is the Einstein tensor and Tab is the "aether 
stress" tensor obtained from varying the aether part of 
the action with respect to the metric. For the metric 
variation we use the inverse metric, but we also have a 
choice whether to consider the independent aether field 
to be a covariant or a contravariant vector. We choose 
contravariant u°' to simplify the stress tensor a bit, since 



the action in the contravariant form has no metric de- 
pendence associated with the C2 and C3 terms in K°'^ran- 
The field equations are thus obtained by requiring that 
the action iQJ be stationary with respect to variations of 
c^'^^ and A. 

The A variation imposes the unit constraint, 

gabu'^u'' = 1. (5) 

The variation gives 

'^aJ°'m - C4UaVmU'^ = Aw™, (6) 

where to compactify the notation we have defined 

= K'^'mn^bu" (7) 

and 

= M°VaU". (8) 

Solving for A using jSJl wc find 

A = u"V, - C4u2. (9) 
The g"'' variation yields the aether stress tensor 

Tab — m{J(a^Ui)-^ — J"^ (a^b) — J(ab)U"^) 

+Cl [(V™Ua)(V"wb) - (V,W™)(VbH™)] 
+CiUaUb 

+ [u„(V„ J"") - C4U^] UaUb 

-^9ab^u- (10) 

In the above expression (O has been used to eliminate 
the term that arises from varying y'—g in the constraint 
term in (Q, and in the fourth line A has been eliminated 
using (PI . The first line contains all of the terms arising 
from varying the metric dependence of the connection. 
Note that it contains terms of second order in derivatives. 
In the last line the notation refers to all of £ in 
except the Ricci scalar term. 

III. SPHERICALLY SYMMETRIC STATIC 
SOLUTIONS 

Our objective in this paper is to consider the weak 
field limit of spherically symmetric static solutions to the 
aether-metric field equations. In spherical symmetry the 
C4 term in the action can be absorbed by the change of 
coefficients 

Cl — > Ci -I- C4 

C3 C3 - C4. (11) 

To see why, note that any spherically symmetric vector 
field is hypersurface orthogonal, hence the twist 

OJa = eabcdu'^^u'' (12) 
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TABLE L Special cases of the action Q. 



Special case parameter values in J^J parameter values in ^ 

General Relativity ci — C2 — C3 — a = t — rj = e = 

equivalent to GR by field redefinition ci+C4 = 0, C2 + C3 = r = e — 

Will-Nordtvedt [1, 4, ^0] C2 = C3 = C4 = rj = e = 

Hellings-Nordvedt \Sj ci + C2 + C3 = 0, C4 = r = 

Einstein-Maxwell-like ci -I- C3 = 0, C2 = C4 = r = rj — 



of the aether vanishes. The identity jl2l | 

= -t^^^w'^+VaWfoVw^- VaMfoVV, (13) 

valid for u satisfying ~ 1, can be used to trade the 
ii^ term in the action Q for an term together with 
the substitution l(TT|l . Since the twist occurs quadrati- 
cally and vanishes in spherical symmetry, that term wiU 
not contribute to the field equations, hence the 0411^ term 
simply modifies the coefficients as indicated in (jllll . Thus 
we henceforth set C4 = without loss of generality, as 
it can be reintroduced at the end via the replacements 
(|ll|l . (Although substitution of the identity (|13(l will not 
change the content of the field equations, it will change 
the value of the lagrange multiplier A for a given solu- 
tion.) 

We have analyzed the asymptotic limit of such solu- 
tions and found at first PPN order the two ERS pa- 
rameters are exactly the same as in pure GR as long 
as C1+C2+C3 0. The special case ci + 02 + 03 = has 
no single characterization. We now describe how these 
results are obtained. 



easily express the field equations in terms of the func- 
tions N{r), B{r), a(r), 6(r) and their derivatives. The end 
result is a set of coupled ordinary differential equations 
coming from both the Einstein equation and the aether 
field equation. Given that there are only three free func- 
tions left after applying the constraint (|16|l . just three 
independent ODE's are needed. We used the aether field 
equation © and the tt and rr components of the met- 
ric equation. The equations are sufhciently complicated 
that it does not seem illuminating to display them here. 



B. Asymptotic weak field limit 

Far from the source, the metric should approach flat 
Minkowski space. In order to examine what happens as r 
approaches infinity we introduce the change of variables 



(17) 



A. Field equations 

A common choice for a weak field analysis is isotropic 
coordinates {t,r,9,(p), which we adopt here. In these 
coordinates the line element is 

ds^ = N{r)dt^ - B{r){dr^ + r^dfl^). (14) 

(Note that r is not the usual Schwarzschild radial coor- 
dinate.) The aether field takes the form 



and the unit constraint becomes 

N{r)a{ry - B{r)h{rf = 1. 



(16) 



The aether field equation © has just t and r compo- 
nents, and the elimination of A reduces this pair to one 
independent equation. 

Solving the field equations is obviously an enormous 
task given the form of the stress tensor so we 

used the symbolic mat h p rogram MAPLE and the Rie- 



tnjp: 

mann tensor package [13. With this package one can 



Around a; = the functions N{x), B{x),b{x) will have 
power series behavior in the form of 



N{x) ^ 1 + Nix + N2X^ + N^x-^ + NiX 
B{x) = 1 + Bix + B2X^ + B3X -I- _D4 
b{x) = &o + bix + b2X^ + h^x^ + 64a;** 



BaX^ 



(18) 
(19) 
(20) 



At this stage it is convenient to use the constraint equa- 
tion H16(l to eliminate a(r) in favor of the radial com- 
ponent b{r). It turns out that asymptotic flatness and 
spherical symmetry generally require the aether to have 
no radial component at infinity (ag = 1, &o = 0) except 
in the Einstein-Maxwell-like case where the action takes 
a special form with an additional symmetry. The first 
order coefficient A'^i determines the Newtonian gravita- 
tional potential, so what we are really interested in are 
the post-Newtonian corrections to this associated with 
the Bi and coefficients. The higher order coefficients 
are post-post Newtonian (and beyond) . Substituting the 
above forms of the functions into the equations of mo- 
tion and performing a series expansion in Maple around 
the point x = ultimately gives a set of algebraic equa- 
tions that can be solved to produce the local power series 
solutions for the fields. 
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C. Series solutions 

To illustrate our methods we first discuss the local 
power series solutions to pure Einstein gravity in isotropic 
coordinates. In this case, the Ci parameters are all set to 
zero and we are left to consider the two coupled ODE's 
for the functions N{x), B(x) given by the vanishing tt 
and rr components of the Einstein tensor. Using the 
procedure described in section IlII Bl we find 



Gt 



-Bi - 2B2 

4 ^ ^ 



iVi ( -Bi - 2B2 



-{Ni + Bi)x^ + 
+ ^Bl - 2N2 - 2B2 



(21) 



(22) 



Thus, at third order in x, the rr equation implies 

Bl = -Ni. (23) 

We can then substitute this result into the fourth order 
equations and solve simultaneously to determine that 



N2 
B2 



1 9 

2^^ 

3 o 
-Nl. 



(24) 
(25) 



For the additional coefficients in B{x) and N{x) we con- 
tinue the process of examining higher order equations, 
substituting in lower order results, and solving simul- 
taneously. These solutions can of course be verified by 
simply expanding the commonly known solution for the 
functions in isotropic coordinates in a power series. 

Now we return to the case of interest and tune the Ci 
parameters back up to non-zero values. At lowest (sec- 
ond) order in x the aether field equation tells us that 



^2(ci+C2+C3)6o(&o + l) =0, 



(26) 



which says that — Q provided that ci + C2 + C3 ^ 
0. This combination of parameters also appears in the 
aether field equation at third order in x, 



-2(ci + C2 + C3)6i 



0. 



(27) 



From equations H26|l and (|27|l it is clear that we have 
two completely different cases depending on whether or 
not ci + C2 + C3 vanishes. These cases must be analyzed 
separately. 



D. Generic case: ci + C2 + C3 7^ 

For this generic case (|26f) shows that asymptotic flat- 
ness of the metric implies that 60 = 0, i.e. the aether 



has no radial velocity at infinity. Together with the con- 
straint this implies that ai — —Ni/2. In addition, H27|) 
tells us that &i = 0. 

Now let us consider the metric equations. The rr equa- 
tion tells us that 



Bl 



(28) 



a result identical to pure GR. We have now determined 
all of the zeroth and first order coefficients in terms of 
iVi, but to examine the higher order ones we must con- 
sider the higher order terms in the expansions of the field 
equations. At fourth order the u field equation is identi- 
cally zero after substituting bi = and Bi = —Ni. Now 
all that remains at this order is to determine B2 and A^2 
using the two Einstein equations at fourth order in x. 
These have the form 



-Nfci -C1N2- 2B2 



-Nl 



InI - 2N2 - 2^2 + ^Nlci 



(29) 
(30) 



Solving these two equations simultaneously yields the fi- 
nal result 



A^2 

B2 



16 



Nlci 



(31) 
(32) 



To determine further coefficients of the power series 
expansion wc move on to consider the field equations at 
fifth and sixth order in x. At fifth order in the u field 
equation we recover the result 



(ci+C2+C3)(62Ari-63)=0 



(33) 



indicating that 62 is a new free parameter in addition to 
A'^i, and 63 = Nib2- The remaining metric equations at 
fifth order are quite complicated so we simply quote the 
final results 



B3 
N3 



16 ^ 96 ^ 

- —Nld 
16 ^ 96 ^ 



(34) 
(35) 



We also examined the sixth order equations to find iV4 
and ,64, but we will not give the results due to their com- 
plexity. However, we note that these coefficients depend 
on both iVi and 62. 

As a final note, we also expanded the equation for 
lambda in ^ and used all of the above results for the 
expansion coefficients to determine at what order lambda 
contributes. This yields 



A = -Nlcix* 
2 ^ 



(36) 



(As mentioned at the beginning of this section, the C4 
dependence of A cannot be obtained via the substitutions 

CD-) 
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E. Special case: ci + C2 + C3 = 

This special case corresponds to the Hehings- 
Nordtvedt theory Q with a unit constraint on the vector 
field. Setting C3 = — C2 — ci from the beginning and re- 
peating the procedure we find that the second order tt 
and rr metric equations imply 

C2bl = 0. (37) 

This special case thus further subdivides into the cases 
C2 7^ and C2 — 0. 



the charged dust interpretation of |9j to the case of ex- 
tremally charged dust.) Thus the exterior solution for 
a star must have A = 0. On the other hand at the ori- 
gin we must have A 7^ to avoid a 1/r singularity in 
the It-field. It does not appear possible to match these, 
so it may be that there are no static spherically sym- 
metric solutions that are regular at the origin. Since the 
Einstein-Maxwell case was already deemed unphysical 
due to the generic appearance of aether shocks, we shall 
not belabor this point here. 

F. Eddington-Robertson-Schiff pEirameters 



1. C2 / 

If C2 is non-zero we again find that u"' has no radial 
component at infinity. There is no single characterization 
of this case. An exceptional sub-case occurs if Ci = 0, 
which falls into the class that is equivalent to GR 
via a field redefinition (and the lagrangian is just R + 
C2RabU°'u^)- In this class the aether field is completely 
unconstrained. If ci 7^ we again find Bi = —Ni as in 
(|^ . while unlike (EH we find 



No 



iiVf(ci-f 6c2-l-4)/(c2 + l). 



(38) 



2. 



C2 







The case C2 = yields the Einstein-Maxwell-like (plus 
C4) sector of the theory, which was previously analyzed 
non-perturbatively in [3. Working through the proce- 
dure for finding the local power series solutions we find 
the Reissner-Nordstrom solution (provided A = 0) with 
bo, 61, and A'^i as free parameters. The freedom appear- 
ing here in bo and bi is a result of the limited gauge 
symmetry 



(39) 



preserving the unit constraint, as discussed in |9|. Specif- 
ically, bi is associated with an "aether charge" while bo 
corresponds to a scaling freedom. This is similar to the 
usual Reissner-Nordstrom case where the general solu- 
tion for the CO- vector potential At is 



Q + Dr 



(40) 



where the D constant is usually set to zero so that the 
field will be at infinity. 

The solutions with A 7^ have the aether aligned with 
the Killing vector, i.e. b{r) = 0. While there always 
exist such solutions in this special case, they are not 
asymptotically fiat except in the even more special case 
ci = — C3 = 2, C2 = 0. (In that case there is a full 
functional freedom in the solution, which corresponds in 



In the usual analysis of the post-Newtonian correc- 
tions to the gravitational field of a static spherical body 
the Schwarzschild line element is rewritten in terms of 
isotropic coordinates and those metric coefficients are 
then expanded to post-Newtonian accuracy. This takes 
the following form for a general gravitational theory : 



1 



2M 



-2(3- 



M 



l~2-f—]\dr^+r^dn^] 
r / 

(41) 

where M is the gravitating mass of the body in geomet- 
ric units and 7 and /? are the Eddington-Robertson-Schiff 
(ERS) parameters of the theory. The parameter 7 mea- 
sures the amount of space curvature produced by a unit 
rest mass and (3 describes the amount of non-linearity in 
the superposition law. 

In the generic case ci -I- C2 -|- C3 7^ 0, we read off from 
equations (|28|l and (|31|l of section IIII Dl that 



7 = 1 
/? = 1, 



(42) 
(43) 



in exact agreement with pure GR. The special case ci -I- 
C2 + C3 = has no single characterization. If ci = it is 
equivalent to GR via a field redefinition, and the aether 
field is arbitrary. If C2 = it is the Einstein-Maxwell- 
like sector, and the exterior is described by the Reissner- 
Nordstrom solution. Hence 7 = 1, but the value of /3 
depends upon the aether charge which is not determined 
by our method. If neither ci nor C2 vanishes then, using 
we find 



7 
(3 



1 

lci_ 
4 



6C2 



C2 



1 



(44) 
(45) 



where was used to obtain /3. This special case cor- 
responds toT = iv = 0^'r] which, as shown by Will 3, 
is dynamically overdetermined in the linearized, uncon- 
strained vector-tensor theory. 



IV. DISCUSSION 

There are two important implications of this analysis. 
First, there appear to be only two free parameters in the 
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local solution around infinity for the generic case C1 + C2 + 
C3 ^ 0, namely Ni and 62- It is possible that analyzing 
the global behavior of the field equations may eliminate 
one of these or demonstrate the existence of even more 
parameters. Based on an analogy from pure GR, the 
metric parameter A^i is determined by the mass of the 
presumed static, central object generating the field. The 
aether parameter 62 cannot be associated with a "charge" 
as in the special case of Einstein-Maxwell due to the 1 /r^ 
fall off. 

The second implication is that in the generic case the 
aether model is quite close observationally to pure GR 
since the ERS parameters match. More precisely, the co- 
efficients of the metric expansions are identical up to B2, 
which differs by a term of relative size ci (or ci -I- C4). 
Other alternative theories of gravity with the same ERS 
parameters are the general vector-tensor theory without 
the unit constraint 14], and the bimetric theories with 
prior geometry of Rosen and of Rastall Q . The fact that 
the ERS parameters are the same suggests that there may 
be a closer relation than might be expected between the 
PPN parameters of unconstrained and constrained vec- 
tor tensor theories. It also is interesting to note the dif- 
ferences between this model and the Brans-Dicke scalar- 
tensor theory. The Brans-Dicke parameters are 



7 
f3 



1 



- OJ 



2 
1 



- OJ 



(46) 
(47) 



where lo is the Dicke coupling constant, which must be 
greater than 500 in order to agree with observation. 

In order to have a comprehensive check on the theory 
in the solar system we need to consider the full post- 
Newtonian approximation scheme. This allows for pre- 
ferred frame effects due to the motion of the solar system 
with respect to the asymptotic preferred frame and is de- 
scribed by ten parameters (two of which are 7 and /3). 
To determine these parameters one perturbatively inte- 
grates the field equations with a fiuid source, imposing 
the condition of regularity at the origin. It seems likely 
that this would fix the value of 62 in the generic case. 
Further tests of preferred frame effects will be found in 
gravitational wave phenomena (briefiy mentioned in 10]) 
such as the orbital decay of binary pulsars, and in strong 
field settings such as black holes. 
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